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ABSTRACT 



Aims. The TeV BL Lac object Markarian 501 is a complex, core dominated radio source, with a one sided, twisting jet on parsec 
scales. In the present work, we attempt to extend our understanding of the source physics to regions of the radio jet which have not 
been accessed before. 

Methods. We present new observations of Mrk 501 at 1.4 and 86 GHz. The 1.4 GHz data were obtained using the Very Large Array 
(VLA) and High Sensitivity Array (HSA) in November 2004, in full polarization, with a final r.m.s. noise of 25 yujy/beam in the HSA 
total intensity image; the 86 GHz observations were performed in October 2005 with the Global Millimeter VLBI Array (GMVA), 
providing an angular resolution as good as 1 10/ias x 40/jas. 

Results. The sensitivity and resolution provided by the HSA make it possible to detect the jet up to ~700 milliarcseconds (corre- 
sponding to a projected linear size of ~500 pc) from its base, while the superior resolution of the 86 GHz GMVA observations probes 
the innermost regions of the jet down to ~ 200 Schwarzschild radii. The brightness temperature at the jet base is in excess of 6 x 10'" 
K. We find evidence of limb brightening on physical scales from S 1 pc to ~ 40 pc. Polarization images and fits to the trend of jet 
width and brightness vs. distance from the core reveal a magnetic field parallel to the jet axis. 

Key words, galaxies: active - galaxies: nuclei - galaxies: jets - BL Lacertae objects: individual: Markarian 501 



1. INTRODUCTION 

The study of extragalactic radio jets is an important area in as- 
trophysics. In radio loud sources, jets contribute a large fraction 
of the total radiated power, and sustain the formation of ener- 
getic kiloparsec scale lobes. While observational properties of 
jets are widely differentiated, they are present in high and low 
power sources, with some common features; on the parsec scale, 
they are relativistic in both types and they are also intrinsically 
iden tical in beamed and misaligned sources. iGiovannini et alj 
(1200 li) have shown that the Lorentz factor in the parsec scale 
jet of low power FRI radio galaxies as well as of more powerful 
FR lis are both in the range F = 3 - 10. With these values, jets 
also appear intrinsically identical in beamed (BL Lac objects) 
and misaligned sources (FRI radio galaxies), if the former have 
j et axis oriented at an average viewing angle of (0) = 18° + 5° 
JGiroletti et al.ll2004bl) . 

In the present paper we focus on the jet structure of the BL 
Lac source Markarian 501. This object is highly active and well- 
studied at all frequencies. Its a ctivity and variabil ity at high en- 
ergy (as high as TeV regime, lOuinn etld] Il996h seems to re- 
quire high Doppler factors and consequently a small angle to 
the line of sight. In the radio band, centimeter VLBI obser- 
vations have revealed a clear limb-brightened structure, begin- 
ning in_the_very_2mierjet, suggestive of a dual velocity struc- 
ture dGiroletti et al.l20()4al hereinafter G04). The complex limb- 
brightened structure makes component identification problem- 



atic and multi-epoch attempts to measure pattern speed con- 
clude that i t is not well defin ed (G04) or in any case at most 
subluminal (Edwards & Finer 2002). These seem to be common 
features in TeV blazars (Finer & Edw ards 2004; Giroletti et aLl 
2006), and theoretical models have been proposed to reconcile 
the m with the ver y high energy emission JGhisellini et al.ll2005l ; 
W ang et al.ll2004l) . 

However, the results obtained so far still leave some major 
questions unanswered. For example, it is not at all clear whether 
the jet velocity structure is intrinsic or produced by the interac- 
tion with the suiTounding medium. We want to understand why 
the properties of the radio jet on parsec scales are different from 
those needed to explain the y-ray emission in Mrk 501, as well 
as in other TeV blazars. Moreover, a change in regime must oc- 
cur on much larger scales, since the large scale structure of the 
source is known to be symmetric rather than one-sided (e.g., 
lUlvestad et al.|[l98 3). We want to investigate if this transition 
is smooth and what is the configuration of the magnetic field in 
the outer jet, which needs sensitive images in total intensity and 
polarization. 

In order to search for an answer to such questions, we need 
to go beyond the capability of the instruments available for ordi- 
nary centimeter wavelength VLBI, which provides information 
only for the region between ~ 1 and ~ 100 milliarcseconds. 
Smaller and larger scale regions remain inaccessible because of 
inadequate resolution and sensitivity, respectively. 
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Fig. 1. Visibility amplitude vs. (m, v)-radius for the GMVA ob- 
servations. 



Improvements in the technical and organizational issues are 
now offering to astronomers VLBI arrays of unprecedented reso- 
lution and sensitivity, such as the High Sensitivity Array (HSAQ), 
and the Global mm- VLBI Array (GMVA0). Thanks to its prox- 
imity and brightness, Mrk 501 is an ideal laboratory for experi- 
ments using these advanced VLBI techniques: it is at z = 0.034 
(1 mas = 0.67 pc, using //q - 70 km s"' Mpc"'); the total flux 
density at 5 GHz is 5 5 = 1 .4 Jy; the Schwarzschild radius for 
its central black hole is estimated around IRs = 10'"^ pc (1.4 x 
lO'^^m as), if we adopt Mbh = 10^ Mq (Rie ger & Mannheiriil 
|2003|) . Using these new facilities, we can therefore access re- 
gions never studied previously: the jet base with the GMVA, and 
the faint, resolved jet region at > 100 mas with HSA. 

In ^ we describe the instruments used for our new obser- 
vations, along with the data reduction methods required by high 
frequency. Results are presented in §3 and discussed in ^ We 
present our conclusions in ^ 



2. OBSERVATIONS 

2.1. High Sensitivity Array observations 

We observed Mrk 501 with the HSA at 1.4 GHz on 26 Nov 
2004. The HSA is obtained by combining in the same array 
the 10 VLBA antennas and other sensitive elements, i.e. the 
Green Bank Telescope (GBT, 100 m.), the phased VLA (27 x 25 
m.), Arecibo (300 m.), and Effelsberg (100 m.). Even without 
Arecibo, whose declination limits do not allow it to observe 
Mrk 501, the collecting area is increased by a factor of 7 over 
the VLBA alone. The sensitivity was also improved thanks to 
a large recording rate (256 Mbps) and a long integration time 
(8hrs). The Effelsberg telescope was in the experiment for the 
first 5 hrs; some failures aff'ected SC, FD, MK, and the VLA 
during part of the observation. 
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Fig. 2. VLA image of Mrk 501. The peak is 1.7 Jy beam ', and 
the contours are traced at (1,5, 10, . . .) X 1.0 mJy beam '. The 
restoring beam is circular with a FWHM of 0.7". Sticks repre- 
sent polarization vectors, with a scale of 1"= 8.3 mJy beam '. 



We reduced the data in the standard way in AIPS, using 
3C345 as a fringe finder, OQ208 as a leakage calibrator, and 
3C286 for the EVR4 calibration. Final images were produced 
both in AIPS and Difmap, with different weighting schemes. 
The source structure is complex, with a strong peak (~ 0.7 Jy) 
and significant diff'use emission (~ 0.8 Jy). Although this pre- 
vents us from reaching the thermal noise, we still achieve in our 
best image a dynamic range as good as 30,000:1, with a noise 
level of ~ 25 fiJy beam"' (lo"). 

During the HSA observations, the VLA (used as phased ar- 
ray) was in the A configuration. We obtained the internal VLA 
data and calibrated and reduced in the standard way to obtain 
also a VLA image of Mrk 501 . 

2.2. Global mm-VLBI observations 

Millimeter VLBI permits a much higher angular resolution 
than ground or space based VLBI at centimeter wavelengths. 
Moreover, it offers the possibility to study emission regions 
which appear self-absorbed at longer wavelengths, with impor- 
tant consequences for our understanding of the physical pro- 
cesses in AGNs in the vicinity of supermassive black holes. After 
years of continuous development and technical improvement, 
the GMVA is now able to provide good quality images in the 
3mm band, with an angular resoluti on of a few tens of micro- 
arcseconds JKrichbaum et alj[2006ah . 

We observed Mrk 501 on 14 Oct 2005 with the Global mm- 
VLBI Array. This experiment tested the sensitivity limits of the 
array, since on the basis of the observed centimeter wavelength 
flux density and spectral index (|G04|) . Mrk 501 was expected to 
be only a few hundred mJy at this frequency. 

The standard frequency was 86.198 GHz, with 16 IFs of 8 
MHz bandwidth each, 2 bit sampling in left circular polarisa- 
tion (LCP). The participating telescopes were Effelsberg, Pico 
Veleta, the Plateau de Bure interferometer, Onsala, Metsahovi, 
and 8 VLBA stations (i.e. all except Saint Croix and Hancock). 
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Fig. 3. High Sensitivity Array image of Mrk 
501. The peak is 762 mjy beam"', and the con- 
tours are traced at (-1, 1,2, 4, . . .) x 0.075 mJy 
beam"' (3cr noise level). As a result of natu- 
ral weighting with a Gaussian taper at 18 MA, 
the restoring beam is 15.1 x 11.4 (mas x mas) 
in PA -11.3°. We mark with letters the region 
of sharp brightness decrease (a) and of limb 
brightening (b) in the extended jet. See Fig. |4] 
for the brightness profile across the jet. 



The European telescopes observed for ~ 9 hours and the 
American ones joined in for the last ~ 6 hours (Mauna Kea only 
for the last ~ 4 hours); the telescopes at Metsahovi and North 
Liberty failed. 

The calibrator 3C345 was readily detected with good signal- 
to-noise ratio. From the fringe fitting of 3C345 we determined 
rates and single-band delays, and applied them to the whole 
data set. We obtained an image of 3C345 and found it to be 
in agreement wit h published images of comparable or slightl y 
lower resolution dLobanov et alj|2000t iLister & Homanll2005l) . 
At this stage, it was then possible to fringe fit Mrk 501 itself, av- 
eraging over the IFs, using a solution interval as long as the scan, 
and setting a SNR threshold of 3.0. Mrk 501 was well detected 
not only between large European apertures but also on baselines 
to the smaller VLBA antennas. Solutions that were obviously 
bad were edited out using SNEDT, and the data were subse- 
quently frequency averaged. Fin al self-calibration and imaging 
were done in Obit (ICottonll2008l) . 

The final amplitude vs. (u, v)-distance plot is shown in 
Fig. [1] The coverage is good in the short baseline range and 
much sparser in the outer part of the (m, v)-plane. Due to the fail- 
ure of the easternmost VLBA antenna, there is also a large gap 
in between the short and long baseline domains. A large base- 
line noise is visible; however, significant emission in the short 
baselines is clearly present. An image of the calibrator 3C345 
and a spectral plot of the resulting phas es vs. spectral channels 
for vis ibilities of Mrk 501 are shown in iGiroletti & Giovanninil 
(l2006h . 



3. RESULTS 

3. 1 . The kpc scale structure 

On kiloparsec scales, Mrk 501 is core dominated with a two 
sided extended structure visible as well, e xtending in PA ~ 45° 
for more than 30" on both sides of the core dUlvestad et al.l 19831 ; 
Kollgaar d et al.lll992l; ICassaro et al.l ll999'). It is straightforward 
to identify this structure with the symmetric extended emission 
characteristic of a radio galaxy and to infer an orientation near 
to the line of sight, in agreement with what is expected from a 
BL Lac source. However, the symmetric emission implies that at 
this distance from the core no relativistic jet remains. 

Thanks to the VLA data available as a byproduct of the HSA 
observations, we obtained a higher resolution VLA image of 
Mrk 501 (see Fig.|2]i. The phased array image is dynamic range 
limited (~ 10000 : 1), and it shows a one-sided emission with a 
short jet like structure in the same PA as the extended symmet- 
ric structure. From this one-sided emission we can derive con- 
straints on the jet velocity (J3c) and orientation (6) with respect 
to the line of sight. At 2" we have yScosfl > 0.36 and at 1", 
/3cos6 > 0.63. This result implies that at 0.67 kpc (projected) 
from the core the jet is still at least mildly relativistic (fi > 0.63). 



3.2. The extended jet 

We obtain a detailed look at the jet of Mrk 501 from the deep 
VLBl observations with the HSA. We show in Fig. [3] a tapered 
image, where baselines longer than 18 MA have been signif- 
icantly down-weighted to increase the signal to noise ratio of 
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Fig. 4. Jet brightness profile across the slice shown in Fig.|3](cen- 
terisatRA 16* 53"' 52^2435035, Dec. +39°45'36".906318,PA 
= -53.4°). 

Table 1. Jet/counter-jet brightness ratio 



the low-surface brightness emission. We achieve a Icr rm.s. of 
~ 25 fiJy beam^' and emission is revealed on the main jet side up 
to a distance of ~ 700 mas from the core, i.e., five times further 
than detected in any previous VLBI observation. No emission is 
detected on the counter-jet side at the level of 3cr noise. In Table 
[T] we give at some selected distances (Col. 1) the jet brightness 
(Col. 2) and the corresponding lower limits to the jet/counter-jet 
brightness ratio (Col. 3). The minimum required velocity (ySmin) 
and largest allowed viewing angles (0max) are then reported in 
Cols. 4 and 5, respectively. The jet opening angle remains con- 
stant {(pj ~ 40°) after the well known bend at 30 mas. 

The jet brightness shows a decrease with increasing jet dis- 
tance from the core; the space distribution is quite uniform, i.e., 
no prominent knot is present in the extended jet. The most note- 
worthy features in Fig. [3] are the relatively sharp brightness de- 
crease at ~ 100 mas (marked with a), which could correspond to 
a shock region, and the jet limb brightening across the slice at 
~ 60 mas from the core (marked by b). The jet brightness profile 
across this slice is shown in Fig. |4] and, although less conspic- 
uous, is similar to the structure visible in the inner jet in higher 
resolution images (G04). 

Besides these features, the jet presents an uniformly dis- 
tributed flux density. If we consider only the region above the 
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Fig. 5. Jet peak brightness and FWHM vs. distance from the core 
obtained from Gaussian fits. A boxcar filter (50 mas) has been 
applied to smooth the data at r > 50 mas. 



3cr noise level, 73.2% of the pixels in the image have a bright- 
ness between 75 and 300 yuJy beam"'. Significant peaks are not 
present and local maxima can be related to small increases in the 
jet emissivity but also to artifacts brought about by the image re- 
construction process. This seems to be a characteristic of the jet 
of Mrk 501 on all scales; even in the inner jet, images with high 
resolution (.G04) tend to show a uniform brightness distribution 
rather than the compact knots observed in other AGN jets. 

In |G04'. we modeled the jet intensity of Mrk 501 as a func- 
tion of jet velocity and radius, us ing the formulas f or an adia- 
batically expanding jet derived bv lBaum et al.l(ll997b in the case 
of relativistic motion. The most sensitive observations available 
in |G04| allowed us to study the jet only to a distance of ~ 100 
mas, still allowing for a degeneracy between magnetic field ori- 
entation and jet velocity. Thanks to the HSA data, we can now 
extend this argument to a distance of almost 500 mas. 

We derived brightness profiles across the jet using the AIPS 
task SLICE on the tapered HSA image for the extended jet, ob- 
taining one slice every 5 mas in PA= -56°. Using the AIPS task 
SLFIT, we fitted single Gaussian components to each profile. We 
show the resulting data as functions of the distance from the core 
in Fig.|5] The fit could be done unambiguously in most cases, al- 
though some slices presented some deviation from a pure single 
Gaussian profile. However, the difference between the area sub- 
tended by the profile and the fit is generally smaller than 5% and 
only a couple of fits (at ~ 200 mas from the core, marked by 
the larger error bars) had to be rejected. At a distance larger than 
~ 450 mas from the core, the best fit FWHM starts decreasing; 
we ascribe this behavior to the insufficient brightness at the jet 
edges, and do not consider any slice at r > 480 mas in our anal- 
ysis. Moreover, we note that for the same reason - and because 
of the limb brightening of the jet - even at smaller distances the 
actual jet FWHM is in some cases larger than the best fit one. 
The implication of this effect are discussed below. In the inner 
jet (r < 30 mas), we keep the fit from G04, whose better resolu- 
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Fig. 6. Mrk 501 at 86 GHz; the restoring beam is 110//as x40yL/as 
in PA -8°. The peak is 45 mJy beam"', and the contours are 
traced at (-1, 1,2,4,8) x4.0mJy beam"'. The Icr noise level is 
~ 1.5 mJy beam"'. The grey scale flux range is -3.0 to 40 mJy 
beam"'. 



tion is important in this region where the jet is not transversely 
resolved in the present image. The best-fit peak brightness and 
jet FWHM have then been smoothed with a boxcar filter (50 mas 
wide, for points at r > 50 mas) to suppress local noise. 

3.3. The core and inner jet structure 

In Fig. |6] we show our Global mm-VLBI Array image of Mrk 
501 at a resolution of 1 10 yuas x 40juas (beam FWHM, PA -8°). 
Mrk 501 is clearly detected at 3 mm and it is dominated by a 
compact, prominent component, ~ 45 mJy beam"' peak bright- 
ness. The visibility data suggest that there is a fair amount of 
extended emission, although the (u, v) coverage is not ideal and 
it is extremely diflicult to image it. 

In our clean image and with modelfitting, we recover ~ 
llOmJy in the core region, including a jet-like feature in PA 
144° and some more diffuse emission in PA ~ -135°. A tenta- 
tive jet knot (~ 7cr) is also visible 0.73 mas south of the core 
(PA 172°). The features in the image plane can be described 
by model-fitting with four Gaussian components. These compo- 
nents are shown overlayed to the 3cr total intensity lowest con- 
tour in Fig. |2l each component is represented by a cross with 
major and minor axis equal to its FWHM, with the major axis 
aligned along the component's position angle. Quantitative re- 
sults from model fitting are reported in Table |2] where r and 6 
are the polar coordinates of the component (re-referenced to the 
core position), bmaj, b^in, and b^ are the deconvolved major and 
minor axis of the component and its position angle, P and / the 
peak brightness and the total flux density. 

Visibility model-fitting in Difmap provides a reduced x^ - 
1.14 with this model. The only quantity that has a significant 
uncertainty in the best-fit (around 10%) is the total flux density 
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Fig. 7. Results from the model-fit with Gaussian components to 
the image plane, overlayed with contours traced at (-4.5,4.5) 
mJy beam"!. The crosses mark the position, major and minor 
axis, and position angle of each component. 



/, while nominal errors on the component positions are typically 
much less than lOjuas. 

Since such values are unrealistically small, we have esti- 
mated independent uncertainties on r taking into account two 
basic parameters for each component: (1) the peak flux den- 
sity P and (2) its compactness. In simple words, the uncertainty 
brought about by noise in the visibility data (i.e. scatter of the 
(m, v)-points) will affect faint diffuse components much more 
than bright compact ones. Therefore, we estimate the uncertainty 
on the position of each component using the following formula: 



Ar = - 

2 



1 V^maj X bn 



P/3o- 



where cr is the image local noise; the formula is therefore 
related to the SNR of the component in such a way that the posi- 
tion of a 3cr feature is not known to better than its mean angular 
radius. The uncertainties reported in Col. (2) are calculated in 
this way and they have been added in quadrature to that on the 
core position (0.004 mas), which is taken as a reference. 

The brightest component, which we identify with the core 
visible at centimeter wavelengths, is still unresolved at 86 GHz. 
We then use our deconvolved size of this component to give an 
upper limit to the dimension of the jet base, and a lower limit 
to its brightness temperature. At z = 0.034, 1 mas = 0.67 pc, 
therefore the deconvolved angular size of the GMVA core cor- 
responds to 0.021 X 0.032 pc. The black hole mass for Mrk 50 1 
is estimated around M^h = 10^ Mq JRieger & Mannheiml2003h . 
which implies a Schwarzschild radius Rs - 1.0 x 10""* parsecs. 
This means that the radio emission originates in a region that is 
smaller than 210 X 320 R5. 

We derive the brightness temperature of this region from the 
following formula: 

B , 

2k 
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Table 2. Deconvolution of Gaussian component fit to the 86 GHz image. 
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(mas) 


(mas) 
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(mjy 
beam"') 


(mJy) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


0.00 


0.0 


0.032 


< 0.048 


169 


40.2 


48.7 ±1.5 


0.08 + 0.03 


-135.6 


0.181 


0.078 


27 


9.7 


51.0 ±4.3 


0.11 ±0.01 


144.4 


< 0.090 


< 0.036 


175 


11.4 


9.0 ±1.1 


0.73 ±0.01 


172.0 


< 0.102 


< 0.030 


170 


12.9 


9.5 ±1.1 



In our observations A - 3.5 mm; moreover, to derive B in 
MKS units, we calculate that 1 beam - 7.52xl0"'^afe ster, where 
a and b are the major and minor semi-axis of the deconvolved 
component, in mas. Therefore: 



r« = 



1.32xlO-i3(fl/7)-'Bnjy 
2x1.38x10-23 



X 12.1 X IQ-^'K 



I.e., 



yiaby^K 



Tb ^5.Sx 10'' xB^jyy 

With the values from Table |2] we find a brightness tempera- 
ture for the core component T^ > 6.8xlO^/ir. If we make the rea- 
sonable assumption that the size of the emitting region is actually 
smaller (e.g. 1/3 of the deconvolved size), we get Tb ^ 6x lO'^A'. 
However, even under this assumption, the result requires a high 
but not extreme value of the Doppler factor at the base of the 
radio jet. In fact, this is not surprising, since the brightness tem- 
perature depends only on the physical length of the maximum 
baseline and on the observed brightness. Therefore, observations 
with similar array at a frequency near the spectral peak can yield 
higher Tb- 

We show in Fig. [8] a spectral plot for both the (average) to- 
tal flux density on kiloparsec scale as measured by single dish 
telescopes, and for the VLBI core, including our new data point 
at 86 GHz. The spectrum of t he V LBI core between 1.6 and 
22 GHz has been presented in IG04| : the core has a turnover at 
about 8 GHz, and then the flux density falls as a power law of 
index a - 0.5. Our new data point follows the optically thin part 
of the core spectrum. Note that as 86 GHz flux density we do not 
just adopt the flux measure of the unresolved core in our image 
but we also include the jet like structure and the diffuse com- 
ponent discussed before; this is because observations at lower 
frequency do not have the angular resolution to distinguish the 
different components. This implies that the turnover frequency at 
~ 8 GHz is related to the whole structure and not to the 86 GHz 
core, whose self-absorption peak is probably located at higher 
frequency. 

We also plot in Fig. [8] (dashed line) the difference between 
the total single dish flux density and the VLBI core one. Apart 
from some fluctuations (due to variability and non simultane- 
ous data), the extended emission from the kiloparsec scale jet 
and lobes region must substantially contribute to the flux den- 
sity even at high frequency. Such extended emission has a rather 
flat spectrum, with index a = 0.3 between 1 .4 and 86 GHz, and 
fluctuations; this odd behavior is likely brought about as a con- 
sequence of variability of the core. 

The detection of a possible jet knot at (r, ff) - 
(0.72 mas, 172°) is in agreement with images at lower frequency. 
Comparing the GMVA data (Table |2]l with modelfits at 15 and 
22 GHz by Edwards & Finer (2002) and G04, the tentative jet 
knot can be identified with the region labelle d as C4, which is 
found at the same distance and position angle. lEdwards & Pined 
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lOwen et alj [19781; iJovce & SimonI II976I: ISteppe et al.l 1 19881) ; 
empty squares show the VLBI core flux: data between 1 .6 and 
22 GHz are mean values taken fromjOOJ; the datum at 86 GHz 
is from the present work. Solid lines connect the points simply 
to guide the eye; the dashed line is the difference between the 
two. Error bars show the standard deviation of averages but for 
the 86 GHz VLBI datum (instrumental calibration uncertainty). 



J2002h found this region to be apparently stationary between 
1995 and 1999, and the positional coincidence lends support to 
this interpretation. At this resolution, the jet has therefore a still 
different orientation with respect to that seen at lower resolutions 
(see G04 and Sect. [lU: the jet PA is ~ 100° at 2 < r < 20 mas 
and ~ 45° at r > 20 mas. 



3.4. Polarization 

In polarized intensity, previous VLBI observations of Mrk 501 
have revealed flux densiti es of a few milliJ ansky, i.e., a few per- 
cent of the total intensity jPushkarev et alJ i2005). Our new HSA 
observations confirm the presence of a significant fraction of po- 
larized flux and reveal interesting details (see Figs. |9] and [TOl i. 
The total flux density in polarization images is almost lOOmJy 
(98.7 mJy), with a peak of 18 mJy beam"' . We plot electric vec- 
tor polarization angle sticks, assuming a Faraday rotation mea- 
sure of Oradm^^. 
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Fig. 9. Polarization images of the inner jet in Mrk 501, overlayed to total intensity contours at 1.4 GHz, traced at (0.2, 0.5, 0.8, 1, 
3, 5, 30, 50, 300) mJy beam"'. Left: EVPA, where 1 mas = 1 mJy beam"'; right: polarized intensity in grey tones between and 5 
mJy beam"'. 



In the inner 10 mas (Fig. |9ll, we have a large cone of po- 
larized flux, with polarization vectors aligned with the jet direc- 
tion. Further downstream, the polarized flux lies predominantly 
toward the southern edge of the jet, in a structure si milar to 
the 'sp ine-sheath' detected in polarization by Pushkare v et al.l 
(l2005l) and in total intensity by G04. At the large bending at 
20 mas, we then find a knot of polarization and a rotation in the 
EVPA position angle. This suggests a strong interaction between 
the radio plasma and an external structure. This interaction could 
be the reason of the change in the source direction, which is fur- 
ther amplified by geometrical eff'ects. 

After the large bend (see Fig.fTOli. the polarization angle be- 
comes well aligned with the jet again at ~ 100 mas from the core. 
At larger distances, the polarized signal becomes weak, with a 
modest preference for a distribution on the south-east side and 
an orientation orthogonal to the jet direction. Finally, the data 
from the VLA only (Fig. |2]i confirm an orientation of the elec- 
tric vector parallel to the jet direction in the extended jet, turning 
slightly clockwise at r ^ 1". 



4. DISCUSSION 

In |3] we have presented our main new results about the core 
and jet of Mrk 501. We now discuss their relevance for our un- 
derstanding of the physics of this source and of AGNs and jets 
in general. 



4.1. The inner core: radio core spectrum and GMVA 
structure 

The nuclear region of Mrk 501 consists of (1) an unresolved 
component: the radio 'core', point-like at our resolution (decon- 
volved size smaller than ~ 30 x 20;uas or 0.020 x 0.014 pc or 
200 X lAQRs), and (2) a faint resolved jet-like structure with a 
large opening angle, similar (taking into account the significant 



difference in flux density an d linear resolution) to the inner struc- 
ture of M87 JLv et al.ll2007l) . 

The unresolved component with a total flux density of about 
45 mJy, is characterized by a relatively low Tb in contrast with 
the higher Tb of M87. However, we note that because of the 
different distance from us of these two sources, our 'core' in- 
cludes the whole M87 s tructure visible in images at 86 GHz 
(iKrichbaum et al.1 l2006bl) . Therefore, the Tb of Mrk 501 is 
the average Tb of a resolved structure where the jet velocity 
could be only mildly relativistic, as suggested b y the detection 
of a counter-jet in the more misaligned M87 JLv et al.l 120071 : 
iKovalev et alj|2007h . 

The diffuse low brightness emission with a total flux density 
of about 65 mJy is interpreted as the continuation of the inner jet 
with spine-shear layer structure. Better sensitivity images are 
necessary to properly map this low brightness feature. 

The lack of a dominant unresolved component is in agree- 
ment with the spectrum shown in Fig. [8] where we can assume 
that the observational data refers to both regions (1) and (2) dis- 
cussed above. Of course, uncertainties related to the variability 
are always present, but because of the regularity of the spec- 
fi-um we can use it to estimate the extension of the radio emit- 
ting r egion, by simply inverting the following formula (Marschen 
Il987h : 



B ^3.2x10-^0* vis: 



?S(l+z)-' 



We assume a local average magneti c field B - .02 G and 
a Dop pler factor i5 = 10 (as discussed bv lTavecchio & Marasch| 
I2OOII) . The low (~ 8 GHz) self-absorption frequency requires 
that the emitting region has a size of the order of 0. 1 mas, in 
agreement with Fig. |6] and that a point-like source, if present, 
is not dominant. Note that the angular size is proportional to 
the magnetic field to the l/4th power, therefore a relatively small 
(even a factor 10) increase of the value of the magnetic field does 
not affect these conclusions. 
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Fig. 10. Polarization vectors overlayed to total intensity contours 
in the extended jet. Contours are traced at (-1, 1,2,4, 8) x 4.0 
mJy beam"', polarization vectors are plotted with a scale of 1.5 
mas mJy \ where unpolarized flux > 2 mJy and polarized flux 
> 0.3 mJy. 



In our GMVA image, we also find a remarkable feature at 
(r, 0) = (0.72 mas, 172°). In fact, lower resolution VLBI obser- 
vations have shown so far that the jet of Mrk 5 01 does not have 
compact jet knots on a few parsecs scale (|G04|) . If the new com- 
ponent is indeed a jet knot, it will be important to re-observe the 
source and test if it has a proper motion that can be followed. 
Alternatively, as suggested by the positional coincidence with a 
feature in previous images, it could be a standing shock at the 
position of a change in jet direction, due to unknown reasons. 
Since this is not far from the first bend in the jet, it may be the 
result of a standing shock from the disturbance causing the bend. 

Moreover, beside this compact feature, a significant amount 
of flux density in the GMVA data remains difficult to image 
and/or modelfit. This is mainly due to the scarce coverage of 
the (m, v)-plane arising from some telescope failures. 

4.2. Jet structure and polarization 

Limb brightening in the jet of Mrk 501 seems to be present on 
scales as small as 0.1 mas, but also after the two main bends at 
~ 2 and ~ 20 mas, where the jet has significantly expanded trans- 
versely. Under a given viewing angle, different Doppler factors 
can arise from diff'erent velocities; therefore, a common explana- 
tion for limb brightening in jets lies in the existence of a velocity 
structure transverse to the jet, with an inner spine and an outer 
shear 

The issue of the presence of transverse velocity struc- 
tures is widely debated at the moment, in the light of both 
analytic al models a nd of numeric simulations (iHarde^ [ 2007t 
Mizuno et al.1 l2007t [Pe rucho et al. 2007; Gopal-Krish na et alj 
2007t iGhisenini et al.ll2OO50 . iChiaberge et all (2000) also sug- 
gested the presence of transverse velocity structures to recon- 
cile observational results with the usually adopted AGN unifi- 
cation model. Direct evidence of limb brightening has recently 



been confirmed, e.g . on the sub-parsec scales in M87 at 43 
GHz (Lv et al."2007') and on a few parsecs scale in 1144-1-35 
(Giovannini et al. 2007). Mrk 501 is unique in the fact that the 
limb brightened structure is visible on scales over three orders of 
magnitude, and in sections of the jet that are diff'erently oriented 
on the plane of the sky. In particular, the limb brightening in the 
mm- VLBI image indicates that a spine/shear structure could be 
present from the very inner part of the jet, close to the region 
where the higher frequency optical emission is produced. Our 
image provides supp ort for the transverse d ual velocity struc- 
ture hypothesized bv IChiaberge et al] (l2000l) on the basis of the 
correlation between the radio and the optical core luminosity in 
radio galaxies and BL Lacs. 

We note that despite the large changes in jet PA, the jet direc- 
tion appears to be constant in time; i.e., we do not have evidence 
of a change in the direction of jet launch, of precession, or of 
variable activity. Only the sharp brightness decrease at ~ 100 
mas from the core could be due to a discon tinuity in the jet ac- 
tivity but the evidence is marginal at most. IConwav & Wrobel 
(1995) tried to interpret this strucutre in terms of a geometrical 
saturating helix model. However, the sensitive HSA image re- 
veals that the emission comes from a conical surface out to more 
that 700 mas, much further out than the transition to cylindri- 
cal surface required in the saturating helix model; moreover, in 
our image as well as in other sensitive low frequency images 
(Giroletti et al. 2004a; Pushkarev et al. 2005), the presence of 
the counter-jet feature described by IConwav & Wrobel (1 19951) 
and predicted by their model is excluded at high confidence 
levels. Clearly, helicity remains an interesting explanation for 
the large position angle changes in the jet ridge line (see, e.j 
the binary black hole system described by IVillata et al.l ll99S 
However, the necessary parameters can not be constrained by the 
available images and the physica l reafity has to be m ore com- 
plex than the geometrical model. iLaing et aP (l2006h have also 
ruled out a helical magnetic field configuration in the nearby, 
low-luminosity radio galaxy 3C 296. The role played by interac- 
tions with the surrounding medium is surely non negligible. 

The polarization structure (Fig.|9]l is also suggestive of spine- 
sheath structures. The most prominent feature in our images is 
the polarized cone departing from the inner region in E-W di- 
rection. We note that because of the lower resolution, this cone 
extends over the full region res olved in the images at higher fre- 
quency ( Pushkarev et al.ll2005h . In our images, polarization vec- 
tors are predominantly aligned with the jet axis, in contrast with 
the dominant p olarization in the external shear perpendicular to 
the jet found by IPushkarev et al.l(l2005h . Because of different fre- 
quency and resolution, a comparison of the datasets is not obvi- 
ous; one can assume that the difference in the polarization vector 
orientation is mainly due to Faraday Rotation or that at 1 .4 GHz 
the dominant polarized flux is from the jet inner spine, and this 
polarized flux has vectors oriented along the jet direction. 

In our image, vectors are aligned with the jet direction also in 
the extended jet structure in NS direction (after the large bending 
at about 50 mas from the core). In this region, the central spine 
is again polarized with vectors aligned with the jet axis and in 
the external shear a marginal evidence of polarized flux oriented 
perpendicular to the jet is present. Only in the bending region 
vectors have a peculiar circular trend but here projection effects 
can be dominant; moreover, in this region a strong interaction 
with the ISM is present and it is the most likely region of the 
change of the jet PA. 

We note that also vectors in the VLA map show a similar 
orientation if we consider that the peak in the VLA polarization 
map is not at the core position but is located in the main pc scale 
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region at about or more 100 mas from the core. We conclude 
that at relatively low resolution and at 1.4 GHz the dominant 
magnetic field structure is perpendicular to the jet axis (E vectors 
aUgned with the jet direction). 

4.3. Jet velocity and orientation 

Our results show that the jet in Mrk 501 is characterized by dif- 
ferent properties on the various scales from a few hundreds to 
several millions Schwarzschild radii. The jet orientation and ve- 
locity, and the ratio between spine and shear contributions must 
significantly change over these scales. It is therefore impossible 
to describe it with constant parameters. 

Since a counter-jet is not detected, however, the jet has to 
be in a relativistic regime even in its faintest and most extended 
region seen in the A configuration VLA observation. The arcsec- 
ond scale structure is symmetric, so we argue that the transition 
to non relativistic velocity has to occur somewhere between pro- 
jected distances of 1 to 10 kpc. 

The jet opening angle remains constant (0j ~ 40°) for sev- 
eral 100 parsecs and any possible re-collimation can take place 
only in the region where the jet becomes confused with the 
noise. Under our estimated viewing angle for this part of the 
jet (0 ~ 15 - 20°), we derive an intrinsic opening angle of 
(f)'. ~ 10- 15°. Jets with larger opening angles are found to have 
lower apparent speeds and Doppler factors in analytical mod- 
els (Gopal-Krishna et al. 2007). This may explain why Mrk 501 
and other TeV bl azars do not ten d to show strong superluminal 
motions (.Finer & Edwardsll2004h . However, it is to be noted that 
the HSA jet of Mrk 501 is several beams wide and does not show 
evidence of any features such as the knot s considered in the an- 
alytical modeling of Gopal-Krishna et al.l (12007) . 

The results from the fits of an adiabatic model to the ob- 
served jet radius and peak brightness presented in Sect. |3.2| can 
be used to constrain the bulk velocity, the orientation, and the 
magnetic field orientation in the various parts of the jet. As we 
have discussed in the section about the polarization properties 
( I4.2l i. we have evidence that the magnetic field is predominantly 
orthogonal to the jet axis. 

First, we recall that froml Giroletti et all (l2004al) it was shown 
that low initial Lorentz factors (F < 5) are ruled out, regardless 
of the magnetic field orientation, since they disagree with both 
the observed limb-brightened structure and the jet/counter-jet ra- 
tio; moreover, they require a jet deceleration between the y-ray 
region and radio jet region that is too strong. 

In Fig. [TT]we use the new HSA data to show the estimated 
jet velocity in the case of parallel and perpendicular magnetic 
field, assuming an initial /3 = 0.998 c and an injection spectral 
index 6 -2{in accordance with lG04l) . In each plot, we draw five 
lines, corresponding to angles to the line of sight of 5°, 10°, ..., 
25° (i.e., in the range of values allowed by the jet sidedness and 
core dominance). 

With these data, we now also rule out models starting with 
F = 15 with magnetic field parallel to the jet axis, since our 
lower limit on the jet/counter-jet brightness ratio near the core 
(R > 4000), and at large distance (R > 140 at 120 mas) is incon- 
sistent with the velocity decrease predicted by a parallel mag- 
netic field adiabatic model at ~ 100 mas. This conclusion is also 
in agreement with the results derived on the polarization prop- 
erties, i.e., a magnetic field in the jet spine orthogonal to the jet 
direction (EVPA parallel to the jet axis). 

The fit with perpendicular magnetic field and an initial 
Lorentz factor F = 15 is in general consistent with the other ob- 



servational constraints. Only in the case of the smallest viewing 
angle (i.e. 6 = 5°) the jet velocity falls ofi^ rapidly after the main 
jet bend; in the extended part of the jet, narrow viewing angles 
are therefore not acceptable. However, it is possible that the jet 
is more closely aligned in its inner part and then it becomes ori- 
ented at a larger 6 after the turn. For all the other viewing angles, 
the fit velocities behave in rather similar ways. After an initial 
decrease (J5 = 0.985 - 0.991, F = 5.8 - 7.5), the jet velocity 
remains relativistic, with small oscillation. This is also in agree- 
ment with the fact that the jet is still one-sided even on scales of 
a few kiloparsecs. 

Finally, we note that the jet FWHM is probably underesti- 
mated at large r. For this reason, we have also tried a fit with an 
input FWHM twice the measured one in the outer jet. The results 
are qualitatively similar to the previous ones: a parallel magnetic 
field is not acceptable and the perpendicular field implies an ini- 
tial deceleration and a more or less constant velocity further out. 
Although the jet velocity is slightly smaller (F = 4 - 6), it still 
remains in the relativistic regime. 



5. CONCLUSIONS 

We have successfully explored new regions in the remarkable 
jet of Mrk 501. Thanks to the great sensitivity of the HSA, we 
reveal that the VLBI jet is one-sided (and therefore in the rel- 
ativistic regime) out to at least 500 parsecs from the core. The 
polarization vectors are clearly aligned with the jet spine, sug- 
gesting that the magnetic field is orthogonal to the jet main axis. 
This is also in agreement with the results of the adiabatic fit to 
the jet brightness and width as a function of distance from the 
core. Limb brightening - already detected on intermediate scales 
by VSOP observations (G04) - is now visible on HSA transverse 
profiles at ~ 60 mas from the core, and is likely present even on 
the sub-parsec scales imaged by the GMVA. 

Despite its presumed weakness, Mrk 501 has in fact been 
clearly detected by the GMVA on sub-milliarcsecond scales. 
This result is encouraging given the performance of the exist- 
ing mm- VLBI array and suggests that not only the brightest 
AGN can be studied on the smallest scales. Present and future 
upgrades of the array (e.g., the installation of new receivers at 
Plateau de Bure, inclusion of new or existing telescopes such as 
the 40 m at Yebes, Spain) are expected to make the instrument 
even more sensitive and reliable. The brightness temperature of 
the most compact component is about 6 x 10'" K. This region 
has a linear size of 0.020 x 0.014pc or, in terms of gravitational 
radii, 200 x lAORs. Significant emission is also revealed in the 
short baselines between the most sensitive telescopes, and awaits 
proper imaging with increased fidelity. 
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